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1.  INTRODUCTION 

The  orotron1'3  is  a  new  type  of  electron  device 
for  generation  and  amplification  of  millimeter 
wave  radiation.  The  device  consists  of  an  electron 
beam  generator  and  collector  and  a  Fabry-Perot 
resonator  containing  one  grooved  mirror  (grating) 
and  one  smooth  mirror.  The  principle  of  operation 
of  the  orotron  is  based  on  the  Smith-Purcell  effect.4 

A  previous  study*  described  the  physical 
principles  of  orotron  operation  and  surveyed  the 
existing  experimental  and  theoretical  literature 
(mostly  from  the  Soviet  Union).  To  augment  the 
information  in  that  report,  quantify  its  conclusions , 
and  aid  in  the  development  of  the  first  operating 
orotron  in  the  United  States  (at  the  Harry 
Diamond  Laboratories  —  HDL),  this  series  of 
reports  will  cover  theoretical  aspects  of  the  orotron 
problem.  A  previous  report'  considered  an 
equivalent  circuit  of  the  orotron  and  derived  the 
output  power  and  frequency  characteristics  of  the 
device.  This  second  report  deals  with  the 
microscopic  mechanism  of  electron  beam 
bunching  in  the  device. 

Earlier  bunching  theories6  7  considered  an 
idealization  of  the  electromagnetic  field  distri¬ 
bution  near  the  diffraction  grating.  This  ap¬ 
proximation  leads  to  an  order-of-magnitude 
estimate  of  the  threshold  current,  but  does  not  give 
the  correct  dependence  on  electron  beam  thickness 
and  grating  parameters.  In  addition,  the  early 
theories  considered  the  motion  of  a  single  electron 
according  to  Newton’s  second  law.  Although  this 
approach  is  correct,  it  is  unwieldy  and  does  not 


?  F  S  Rusm  and  G  O  Bogomolov.  Soviet  Patent  / 95557  ( 1965). 

*  ,V  Tagushi  and  S  Ono.  Report  of  Research  Group  on  Electron  Devices. 
Sendai.  Japan  (March  1964) 

1  I  P  Shestopalov.  Diffraction  Electronics.  Khar'kov.  Moscow  (1976) 
(Trans  hv  L  S.  Joint  Publications  Service.  April  I97g). 

4  S  J  Smith  and  E.  M  Purcell,  Phys  Rev.  92  (1953).  1069 

'  R  P  Leavitt.  Characteristics  of  Orotron  Oscillation  and  Amplification.  I 
Power  and  Frequency  Characteristics,  Harry  Diamond  Laboratories  HDl.-TR 
1X99  (July  1979) 

h  M  B.  Tsevthn.  G.  A  Bemashevskiy.  V  D  Kotov,  and  I  T  Tsitson .  Radio 
Eng  Elect  Phys  .  22  (1977),  132 

7  K  Mizuno.  S.  Ono.  and  Y  Shibata.  IEEE  Trans  Electron  Dev..  F.D-20 
0973).  749 

*  D  E  Wortman  and  R  P  Leavitt.  Near  Millimeter  Wave  Orotron  Research 
Study.  Ham ■  Diamond  Laboratories  (draft) 


lend  itself  to  straightforward  generalization  to  a 
nonlinear  theory. 

We  consider  the  bunching  problem  from  the 
point  of  view  of  distribution  functions  that  describe 
the  electron  beam  as  a  whole,  and  we  consider  the 
time  evolution  of  the  beam  via  the  Boltzmann 
equation.  This  method  provides  a  connection 
between  the  orotron  and  other  free-electron  lasers 
that  have  been  analyzed  by  this  method.89  In 
section  2,  a  general  formula  is  derived  for  bunching 
of  an  electron  beam  by  an  applied  longitudinal  rf 
electric  field.  Section  3  applies  this  result  to  the 
orotron  and  computes  the  output  power  of  the 
device;  by  equating  the  power  radiated  by  the 
electron  beam  to  the  power  absorbed  by  the  open 
resonator,  we  derive  an  expression  for  the  starting 
current  for  the  device.  Section  4  justifies  the 
representation  of  the  electron  beam  as  a  complex 
frequency-dependent  admittance  (in  the  language 
of  our  earlier  report5)  and  derives  the  electronic 
tuning  transconductance  of  the  orotron. 

2.  ELECTRON  MOTION  IN 
LONGITUDINAL  ELECTRIC 
FIELD 

In  the  orotron  ( as  opposed  to  other  free-electron 
laser  devices),89  electrons  are  bunched  by  the 
longitudinal  rf  electric  field  (the  component  of  the 
electric  field  along  the  direction  of  electron 
motion)  in  the  open  resonator.  In  this  section,  we 
analyze  the  motion  of  an  electron  beam  in  a 
monochromatic  longitudinal  electric  field,  subject 
to  the  following  assumptions: 

a.  The  beam  is  immersed  in  a  longitudinal  dc 
magnetic  field  that  is  intense  enough  so  that 
electron  motion  transverse  to  the  beam 
propagation  direction  can  be  neglected.  This 
assumption  allows  us  to  consider  electron 
motion  along  one  direction  only,  which  we 
call  the  y-direction. 

'  R.  P  Leavitt.  Characteristics  of  Orotron  Oscillation  and  Amplification.  I 
Power  and  Frequency  Characteristics.  Ham  Diamond  Laboratories  HDL-TR 
1X99  (July  19  79) 

*  F.  A  Hop/  P  Mevstre.  G  T  Moore,  and  M  O  Scully,  in  Physics  of  Quantum 
Electronics.  Yol  5.  Novel  Sources  of  Coherent  Rodiation.  Addtson-WesU 
Publishing  Co .  Inc.  Reading.  MA  (I97&) 

4  P  Sprang le  and  A  Drobot.  IEEE  Trans  Microwave  Theory  Tech.,  MTT-25 
(1977).  52* 
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b.  The  electron  beam  accelerating  voltage  is  time  t  is 
low  enough  so  that  relativistic  effects  can  be 

neglected.  This  assumption  allows  us  to  use  «, 

the  nonrelativistic  equations  of  motion  to  p(y,t)  =  — e  (  dv  f(y,v,t)  ,  (1) 

analyze  the  electron  beam  propagation.  ^  —  «. 


c.  The  density  of  electrons  in  the  beam  is  low 
enough  so  that  the  forces  of  repulsion  be¬ 
tween  the  electrons  can  be  neglected.  This 
assumption  avoids  the  need  to  use  a  self- 
consistent  solution  of  Maxwell's  equations 
and  the  dynamical  equations. 

d.  The  longitudinal  electric  field  in  the 
resonator  is  weak  enough  so  that  it  may  be 
considered  as  a  perturbation  on  the  freely 
propagating  electron  beam.  This  assump¬ 
tion  allows  us  to  make  a  perturbation 
expansion  of  the  dynamical  equations  and  is 
consistent  with  the  point  of  view  that  regards 
the  orotron  as  a  “free-electron  laser  of  the 
Smith-Purcell  type.”8 

One  way  to  approach  the  problem  of  bunching 
in  orotrons  is  to  consider  the  behavior  of  a  single 
electron  under  the  influence  of  the  longitudinal  rf 
field  of  the  open  resonator.67  This  approach  is 
cumbersome  and  leads  to  the  same  results  as  those 
given  here.  In  this  work,  we  approach  the  problem 
via  the  distribution  function  approach,  which  is 
equivalent  to  the  previous  method.  Furthermore, 
the  distribution  function  approach  can  deal  with 
increased  density  systems  where  the  single  electron 
approach  would  be  hopelessly  complicated. 

The  distribution  function  describes  the  density 
of  electrons  in  coordinate  and  velocity  space.  To 
be  more  precise,  if  f(y,v,t)  is  the  distribution 
function,  then  fiy,v,t)  dy  dv  is  the  number  of 
electrons  having  coordinates  between  y  and 
y  +  dy  and  velocities  between  v  and  v  +  dv  at 
time  t.  In  terms  of  the  distribution  function, /he 
charge  density  of  the  electron  beam  at  position  y  at 

6M.  B.  Tseytlin.  G.  A.  Bernashevskiy,  V  D  Kotov,  and!  T.  Tsitson',  Radio 
Eng  Elect  Phys.  22  (1977).  132. 

7  K.  Muuno.  S.  Ono,  and  Y.  Shtbata.  IEEE  Trans.  Electron  Dev..  ED-20 
0973).  749. 

8  F.  A.  Hopf.  P.  Mevstre.  G  T  Moore.  andM.  O.  Scully,  in  Physics  of  Quantum 
Electronics.  Yol.  5.  Novel  Sources  of  Coherent  Radiation.  Addison-  Wesley 
Publishing  Co..  Inc..  Reading.  MA  (1978). 


and  the  current  density  at  position  y  at  time  t  is 


Joo 

dv  v  f(y,v,t)  .  (2) 

-  CO 

In  equations  ( 1 )  and  (2),  e  is  the  magnitude  of  the 
electron  charge.  The  dynamical  equation  of  mo¬ 
tion  for  the  electron  beam  in  terms  of  its  distribu¬ 
tion  function  is  given  by  the  Boltzmann  equation:10 


df  _df _ e  p  df 

dt  V  dy  m  dv 


(3) 


where  m  is  the  electron  mass  and  E  =  E(y,t)  is 
the  longitudinal  rf  electric  field.  By  integrating 
equation  (3)  over  velocities  and  using  equations 
( 1 )  and  (2),  the  reader  may  verify  that  the  charge 
conservation  equation 


A 

dy 


,  dp  = 

dt 


is  satisfied. 

In  deriving  equation  (3),  we  have  used  assump¬ 
tions  a  to  c.  That  is, 

a.  Equation  (3)  is  one-dimensional. 

b.  Equation  (3)  is  nonrelativistic. 

c.  Equation  (3)  does  not  contain  the  electronic 
repulsion  forces. 

Within  these  constraints,  equation  (3)  exactly 
represents  the  electron  motion. 

Assumption  d  is  incorporated  into  our  discus¬ 
sion  by  expanding  equation  ( 3)  in  a  power  series  in 
the  strength  of  the  electric  field.  We  introduce  the 
strength  parameter  X  and  let  E  —  XE  in  equation 
(3).  Then  we  expand  the  distribution  function  in  a 


R.  C.  Davidson.  Methods  in  Nonlinear  Plasma  Theory  .  Academic  Press. 
Sen  York  (1972) 
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power  series  in  A: 


fl(y,v,t)  =  2  A’f,,  (y,v,t)  (4) 

n  =  0 

This  expansion  is  then  substituted  into  equation 
(3),  and  the  coefficient  of  A"  is  equated  to  zero.  The 
result  is,  for  n  =  0, 


df0  df0 
~dt  +  y  W 


=  0 


and,  for  n  >  0, 


(5) 


,  dfn  _  _C_  p  dfn  —  1 

dt  v  dy  m  dv 


(6) 


We  choose  the  zeroth-order  distribution  func¬ 
tion  to  be  a  function  of  velocity  only,  to  correspond 
to  the  physical  situation  of  a  uniform  dc  electron 
beam.  The  normalization  of  f0  is  chosen  such  that 
the  zeroth-order  current  density  is 

jo  =  -e  C  dv  V  f0  (v)  =  Io  /A0  ,  (7) 

J  —  oo 


where  Iq  is  the  total  beam  current  and  Ao  is  the 
beam  cross-sectional  area.  This  choice  of  the 
zeroth-order  distribution  function  satisfies  equa¬ 
tion  (5)  since  f0  is  independent  of  y  and  t. 


The  longitudinal  rf  electric  field  is  chosen  to  be 
monochromatic  and  can  be  written  in  the  form 


E(y,t)  =  E0(y)e~'.<ul  +  c.c.  ,  (8) 

where  w  is  the  angular  frequency  of  the  field  and 
c.c.  means  complex  conjugate.  We  wish  to  con¬ 
sider  the  first-order  contribution  to  the  distribution 
function;  equation  (6)  becomes,  for  n  =  1, 


Equation  (9)  is  a  linear,  first-order,  ordinary 
differential  equation,  which  may  be  solved  by 
standard  methods."  The  result  is 

«.(V.v)  =  ~  eiujy/v  fo»  X 


Cy  dy'  e-'^Eofy')  ,  (11) 

*)  —  oo 

where  we  have  used  the  condition  a,  ( —  °°  ,v)  =  0. 
Substituting  equation  (11)  into  equation  ( 1 0)  gives 
the  first-order  distribution  function. 


The  first-order  contribution  to  the  current  den¬ 
sity  is  obtained  from  equation  (2).  The  velocity 
integral  in  equation  (2)  may  be  integrated  by  parts, 
giving  the  result 


ji(y)  =  - 


ie2io 

m 


f  ei  toy/v  f(  (v)  x 

J-«,v 


(12) 

dy'  (y  -  y')e~iuty/v  E0(y')e~"'"  +  c.c. 


The  power  radiated  by  the  electron  beam  is  given 
by 


P  rad 


I 


jE  dr  ^ 


(13) 


where  the  integral  is  over  the  volume  of  the 
electron  beam  and  the  brackets  indicate  a  time 
average.  Substituting  equations  (8)  and  (12)  into 
equation  ( 1 3)  yields  the  result 


Prod 


e2w 

m 


X 


J”  beam  * 


dy^_oo  dy'  X 


v  “dy1  ~  iw0fl  =  "m  ’  W 

where  f^(v)  =  df0/dv.  In  writing  equation  (9), 
we  have  introduced  the  first-order  distribution 
amplitude  a,  via 


tiw(y  -  y’)/v  Eo(y)E„(y')  +  c.c. 


(14) 


where  the  integral  over  the  transverse  coordinates 
x  and  z  is  over  the  cross-sectional  area  of  the  beam. 
The  longitudinal  field,  Ep(y ) ,  depends  implicitly  on 
x  and  z  in  addition  to  y. 


fi(y,v,t)  =  o,(y,v)e  “■*  +  c.c. 


(10) 


11 H  W  Reddick ,  Differential  Equations ,  John  Wilev  and  Sans,  Inc  ,  Nf* 
York  ( 1 94  J>. 
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Equation  (14)  allows  the  computation  of  the 
power  radiated  by  the  electron  beam  provided  that 
the  distribution  of  electric  field  strength,  E0(y),  is 
known.  It  is  more  convenient,  however,  to  charac¬ 
terize  the  electric  field  in  terms  of  its  Fourier  trans¬ 
form: 


This  is  the  main  result  of  this  section;  it  allows  the 
computation  of  the  power  radiated  by  an  electron 
beam  immersed  in  a  longitudinal  monochromatic 
electric  field.  In  the  next  section,  we  consider  the 
orotron,  in  which  the  electric  field  is  quasi- 
periodic. 


Joo 

_oo  e  ~ipy  E0(y)  dy  .  (IS) 

In  terms  of  the  Fourier  transform,  the  double 
integral  in  equation  (14)  becomes 


f  *  r> 


ei<v-/yvE.(y)Eo(y-)  = 


2 


2 


dp  |e0(p)|  2 


to 


,  (16) 


where  P  means  principal  value.  By  using  equation 
(16),  equation  (14)  becomes 


In  what  follows,  we  are  interested  in  the  radia¬ 
tion  from  an  electron  beam  in  which  all  of  the  elec¬ 
trons  are  traveling  at  a  velocity  v0.  By  using  the 
normalization  given  by  equation  (8),  the  distribu¬ 
tion  function  is 

f,(v)  =  -  — S  (v  -  v„)  .  (18) 

?vr>Aj> 

and  the  power  radiated  at  frequency  w  is 


Prod 


ewl0 

mvoA„ 


X 


l  dx  dz  - — 

(beam  aw 


■  09) 


3.  RADIATION  BY  ELECTRON 
BEAM  NEAR  PERIODIC 
STRUCTURE 

In  applying  equation  (19)  to  radiation  in  an 
orotron,  we  must  account  for  the  detailed  structure 
of  the  electromagnetic  field  distribution.  Consider 
the  schematic  diagram  of  the  orotron  shown  in 
figure  1.  The  electron  beam  passes  close  to  the 
reflecting  diffraction  grating,  and  the  interaction  of 
the  beam  with  the  longitudinal  rf  field  of  the 
resonator  (as  modified  by  the  grating)  causes  the 
beam  to  radiate,  as  shown  in  section  2. 


The  electromagnetic  field  of  the  open  resonator 
of  the  orotron  may  be  considered  as  a  natural  mode 
of  the  resonator,  with  a  fine  structure  imposed  by 
the  grating.  This  fine  structure  causes  the  radiation 
characteristic  of  the  orotron. 

Consider  a  plane  wave  of  frequency  w  incident 
normally  on  the  surface  of  the  grating,  as  shown  in 
figure  2.  Let  the  period  of  the  grating  be  f .  The 
electromagnetic  field  may  be  characterized  by  the 
x-component  of  magnetic  induction,  given  above 
the  grating  by* 


Hr  = 


2  cos  k  (z  —  z0)  +  £  a,  X 


r  =  1 


cos 


M 


e  'r ‘ 


+  c.c. 


(20) 


where  k  =  w/c,  c  is  the  speed  of  light,  H0  is  the 
amplitude  of  the  incoming  wave,  and 


rr  = 


(¥)'- 


(21) 


*D.  H  Wortman  and  R  P.  Leavitt,  Near  Millimeter  Wove  Orotron  Research 
Study,  Harry  Diamond  Laboratories  (draft). 
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SMOOTH  REFLECTING  MIRROR 


ELECTRON  COLLECTOR 


ELECTRON  BEAM 


9 


E0(y)  =  - 


1H0 

ce„ 


Equation  (20)  represents  the  most  general  free- 
space  solution  of  Maxwell’s  equations  corre¬ 
sponding  to  an  incoming  plane  wave  that  satisfies 
the  periodicity  requirement 


2  sin  k  (z  —  z0)  + 


Hx  (y  +  l.z)  =  Hx  (y,z)  (22) 


imposed  by  the  periodic  grating.  The  unknown 
coefficients,  arin  equation  (20),  are  determined  by 
solving  the  appropriate  boundary  value  problem  at 
the  grating  surface.  We  do  not  concern  ourselves 
with  the  solution  of  the  boundary  value  problem 
here,  but  we  consider  the  coefficients  a,  as  being 
known.*  The  y-component  of  the  electric  field  is, 
from  Maxwell’s  equations, 

=  J_  = 

‘  O)£0  dz 


iHo  ' 

C£„ 

COS 


1  f 

2  sin  k  (z  —  z0)  +  ar  —  X 

'  r  =  \  k 


/  2r7ry  \  . 

rr)e  r 


+  c.c.  ,(23) 


X  fx(x)f,.(y)  .  (25) 

This  is  the  electric  field  that  is  used  in  equation 
( 1 9)  to  compute  the  radiated  power. 

In  computing  the  Fourier  transform  of  this  elec¬ 
tric  field,  it  is  necessary  to  make  a  few  approxima¬ 
tions  whose  validity  depends  on  the  fact  that  the 
envelope  function  f,.(y)  is  slowly  varying  com¬ 
pared  with  the  variation  of  the  terms  in  brackets  in 
equation  (25).  The  function  varies  slowly  if  the  di¬ 
mension  of  the  resonator  mode  alongy  is  very  large 
in  comparison  with  the  period  of  the  grating.  By 
using  equation  (15),  the  Fourier  transform  be¬ 
comes 


where  e0  is  the  vacuum  permittivity. 


E0(P)  =  - 


iH„  ~ 

C£«  N  =  - 


I 


1/2 

-1/2 


dy  X 


The  electric  field  of  equation  (23)  corresponds 
to  an  incident  plane  wave;  it  must  be  modified  to 
account  for  the  mode  pattern  in  the  open  resonator. 
For  a  large  class  of  open  resonators,  the  mode 
pattern  assumes  the  form  fx(x)fv(y),  where  f*  and  fv 
are  Gaussians  or  Gaussians  multiplied  by  Hermite 
polynomials  describing  the  variation  of  electro¬ 
magnetic  fields  within  the  resonator.  The  functions 
are  normalized  to  unity: 


e  ~ip(y  +  wfx(x)f,(y  +  Nl)X 


2  sin  k  (z  —  zcl)  4- 


2nTy 


.  (26) 


_  .  dx  |  fr(x)  |  2  =  1  ,  (24a) 

dy  |  My) |  2  =  1  ■  (24b) 

In  terms  of  these,  the  longitudinal  electric  field  of 
equation  (8)  becomes 

The  determination  at  the  coefficients  a,  will  be  discussed  in  a  forthcoming 
report  in  this  ser-es 


where  we  have  broken  the  y-axis  up  into  segments 
of  length  l  and  used  the  periodicity  condition, 
equation  (22). 

It  is  known  that  radiation  from  the  grating 
should  peak  at  frequencies  near  those  given  by  the 
Smith-Purcell  condition4  * 


4  S  J  Smith  and  ft*.  M  Purcell.  Phys  Rev.  92  (I9SJ).  1069 
*  D  E.  Wonman  and  R  P  Leavitt.  Sear  Millimeter  Ware  Orotwn  Research 
Study.  Horn-  Diamond  Laboratories  (draft). 
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-  -  —  " 


Therefore,  referring  to  equation  (19),  we  shall 
consider  the  Fourier  transform  near  p  ~  2rrr/t, 
where  r  is  an  integer.  *  Assuming  that  fy  does  not 
vary  appreciably  over  an  interval  of  length  1,  we 
have 

~  / 2nr  + 


fx(x)  2  fy(N t )  t~iSN  X 

C£0  N  =  —  <• 

cm  r 

I  e  2  sin  k  (z  -  z0)  +  (28) 

“  I>  /2r'rry\  -r,/ 
i  a/  — —  cos  l  —  I  e 

/  =  l  k  V  1  '  -1 

The  integral  selects  out  the  term  with  r'  ~  r;  the 
remaining  sum  can  be  approximated  by  an  integral. 
The  result  is 


~  2nr  +  8 

En  I - T— 


^  *  6  X 
JL  {(s)2 

d8  >\l} 


at  a  frequency  i o  =  (v0/l)  (2rrr  +  8). 

Specifically,  equation  (31)  may  be  evaluated 
for  the  case  of  a  TEM20  mode  in  the  resonator.  *  In 
this  case,  we  have 


(r(x)  = 


4  (2j r) 1,2  w; 


H,  e  ~x2/»i 


where  H2  is  the  second-degree  Hermite  poly¬ 
nomial,  and 


fv(y)  = 


'■v2/wy  .  (33) 


^  5  (7 


,  (29) 


The  width  parameter,  w*.  is  chosen  equal  to  the 
width  of  the  electron  beam,  D.  The  integral  over 
the  electron  beam  cross  section  in  equation  ( 3 1 )  is 
then 


1  dx  dz  fT(x)  e 

/beam  1 


f>)  =  J 


dy  e  ~ipy  fTy)  (30) 


is  the  Fourier  transform  of  the  envelope  function, 
fv.  We  find  the  radiated  power  to  be,  from  equation 
(19), 

=  _  ell0  |H0l2a2r2 
rad  4mv^wA0eo 


*  We  refer  to  an  orotron  radiating  at  frequencies  near  2rv.  v()  ft  or  an  orotron 
operating  on  the  nh  harmonic 


,  (34) 


wfiere  a  is  the  height  of  the  beam  above  the  grating 
and  A  is  the  beam  thickness  (A0  =  DA).  The 
Fourier  transform  of  f,.(y)  can  be  evaluated  in  a 
straightforward  manner  to  give 


~  ( 8\  ,  u  i ' 

f,.  ( j  )  ^  (2^  (w,)*  e  '  .  (35) 

The  quantity  |  H0 1 2  in  equation  (31)  may  be 
eliminated  in  terms  of  the  total  energy  in  the  open 
resonator, 

*  71*  TF.M 20  node  is  the  mode  most  common! \  encountered  tn  practice 


Thus,  the  starting  current  is  given  by 


U  =  ~  <  J(eoE2  +  MoH2)  dr>  ,(36) 

where  the  integral  is  over  the  volume  of  the  open 
resonator  and  juo  is  the  vacuum  permeability. 
Neglecting  the  contribution  of  the  evanescent 
waves  in  equation  (25),  we  obtain 


Ij  (min) 


mvJoj2DL£0 
0.0095eaJ  f2c2Wy  Q 


U  = 


4LHp 

£0C2 


(37) 


where  L  is  the  mirror  spacing  in  the  open  resonator 
and  c  is  the  speed  of  light  in  vacuo.  Inserting 
equations  (34),  (35),  and  (37)  into  equation  (31) 
gives  the  result 

=  _  Q.2(27r)l/2eI0Ua2rV 
rad  16  mv2(uDL£o 


Equation  (43)  has  been  evaluated  for  the 
orotron  being  built  at  HDL  to  operate  in  the 
fundamental  mode  (r  =  1)  at  75  gigahertz.  The 
electromagnetic  boundary  problem  in  equation 
(20)  has  been  solved*  with  t  =  0.4  millimeter 
to  give 


a,  =  0.2796, 

r,  =  1.563  X  1  OVmeter.  (44) 


e 


-2Vra 


2Tr  A 


8 


-W2VS2/ 2P 
e  y 


•(38) 


We  consider  the  case  of  a  beam  that  grazes  the 
grating;  therefore,  a  =  0.  The  other  parameters  of 
the  system  are 


The  orotron  actually  radiates  at  a  frequency  given 
by  the  peak  of  equation  (38);  the  peak  is  at 


Thus,  the  peak  output  power  from  the  device  is 


Pra</(max)  =  0.0095  |lo|  U 


ea2T?c2w2 

mvjJwDLe,, 


X 


The  starting  current  for  orotrons  may  be  deter¬ 
mined  by  requiring  that  the  power  gain  given  by 
equation  (41)  equal  the  power  loss  in  the  resonator, 
given  by 


^  loss 


Uoj 

~Q 


(42) 


m  =  9.1  X  10  ■'  kilogram, 
v0  =  3  X  107  meters/second, 
uj  =  2tt  X  75  X  109/ second, 
a  -  0, 

A  =  0.03  centimeter, 

D  =  1  centimeter, 

L  =  2  centimeters, 

£0  =  8.854  X  10“ 12  farad/meter, 
e  =  1.6  X  10  19  coulomb, 
c  =  3  X  108  meters/second, 

Wj,  -  1  centimeter, 

Q  =  5000  (estimated).  (45) 

*  The  determination  of  the  coefficients  ar  will  be  discussed  in  a  forthcoming 
report  in  this  series 
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Using  these  parameters,  we  obtain 

Ij  =  34.5  milliamperes.  (46) 


ip  r  dpi  e0(p)1  ; 
;  J"  '-£) 


4.  ELECTRONIC  TUNING  CHARAC¬ 
TERISTICS  OF  OROTRON 

In  a  previous  report  of  this  series,5  we  showed 
that  the  frequency  characteristics  of  an  orotron 
could  be  derived  from  an  equivalent  circuit 
representation.  The  electron  beam  was  repre¬ 
sented  in  this  circuit  by  a  negative  conductance 
that,  if  sufficiently  large  in  magnitude,  acted  as  a 
feedback  element  leading  to  oscillation  of  the 
circuit.  This  model  predicted  the  correct  behavior 
of  the  orotron  with  regard  to  the  threshold 
condition  and  the  behavior  of  the  output  line  width 
below  and  above  threshold.  One  of  the  parameters 
in  the  model,  the  threshold  current,  Ij,  has  been 
calculated  in  this  report  (sect.  3)  from  first 
principles.  However,  the  electron  beam  is  not 
actually  represented  by  a  frequency-independent 
conductance,  as  was  assumed  in  the  previous 
report. 


where  fi  is  a  proportionality  constant  and  we  have 
considered  the  sharply  peaked  distribution  func¬ 
tion  given  by  equation  (18).  We  have  chosen  the 
sign  in  equation  (47)  so  that  real  (Re)  Y  beam  >  0 
corresponds  to  radiation. 

Section  3  is  concerned  with  the  evaluation  of 
the  first  term  in  equation  (47).  Since  at  this  point 
we  are  interested  only  in  the  frequency  depen¬ 
dence,  we  can  represent  the  result  of  that  calcula¬ 
tion  by  a  “normalized  conductance”: 

Re  ybeam  =  -  u  e<'  '  “J>/2  .(48) 

Ii(min) 

where 

wv5  /  to  2rrr 


Consider  equation  ( 1 4)  for  the  radiated  power. 
The  real  part  of  equation  (16)  represents  the  effect 
of  that  part  of  the  current  that  is  in  phase  with  the 
applied  electric  field;  this  is  the  only  part  that 
contributes  to  the  radiated  power.  However,  the 
imaginary  part  (which  does  not  contribute  to  the 
time-averaged  power)  representing  the  out-of¬ 
phase  (by  90  deg)  current  component  can  be 
considered  as  being  due  to  an  “imaginary 
conductance”  (rigorously,  a  susceptance). 
Therefore,  the  overall  effect  of  the  electron  beam  is 
to  introduce  a  frequency-dependent  complex 
admittance  (rather  than  a  simple  conductance,  as 
was  assumed  earlier)5  into  the  equivalent  circuit. 


This  result  is  obtained  by  squaring  equation 
(35)  and  differentiating  with  respect  to  S.  The  nor¬ 
malization  is  chosen  so  that,  at  u  =  —  1  (where 
Re  y  beam  attains  its  maximum)  and  at  a  current 
10  =  15  (min),  the  beam  conductance 
Re  y  beam  ~~  I • 

From  equation  (47),  the  imaginary  (Im)  part  of 
the  beam  admittance  is  given  by 

Im  y beam  (u)  = 

du  Re  y  beam  (u) 


We  consider  an  admittance  that  is  proportional 
to  the  derivative  of  equation  (16): 


This  integral  may  be  performed  easily  by  using 
equation  (48);  the  result  is 


*  R.  P  Leavitt.  Characteristics  of  Orotron  Oscillation  and  Amplification.  I 
Power  and  Frequency  Characteristics.  Harry  Diamond  Laboratories  HOI.  Tfl 
1899  (July  1979) 


Im  y  beam  (n)  = 
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where  F  is  a  confluent  hypergeometric  function. 
The  real  and  imaginary  parts  of  y  beam  are  shown  in 
figure  3  in  the  region  u  <  0,  where  the  beam 
radiates. 


y  =  ^1  -  Re  ybeam  (u)J  - 

[to  —  w0  "1 

2Q  - +  Im  y^m  (u)J  ,(52) 


Consider  the  orotron  equivalent  circuit  shown 
in  figure  4,  which  is  the  circuit  given  earlier5 
generalized  in  accordance  with  the  above  in  terms 
of  normalized  admittances.  The  overall  admit¬ 
tance  of  the  circuit  (normalized)  is  given  by 

y  Yo  +  Sour  —  Ybeam  *  (51) 

where  y0  is  the  normalized  admittance  of  the  open 
resonator  and  goul  is  the  normalized  conductance 
representing  the  output  coupling.  By  using  the 
results  of  reference  5,  equation  (5 1 )  may  be  written 
as 


where  Q  is  the  quality  factor  of  the  passive  open 
resonator  and  w0  is  the  resonant  frequency  of  the 
resonator.  The  normalization  is  chosen  so  that 
Re  y0  +  gou(  =  1. 

If  we  consider  the  orotron  as  an  ideal  oscillator 
(that  is,  we  neglect  the  noise  source  in  fig.  4).  the 
condition  for  oscillation  is  that  both  the  real  and  the 
imaginary  parts  of  equation  (52)  must  vanish.  That 
is, 

Re  y beam  =  1  (53) 

for  the  real  part  and 


5  R.  P.  Leavitt.  Characteristics  of  Orotron  Oscillation  and  Amplification.  /. 
Power  and  Frequency  Characteristics.  Harry  Diamond  Laboratories  HDL  TR 
1899  (July  1979). 
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Figure  4.  Equivalent  circuit  of  orotron. 


for  the  imaginary  part.  By  using  equation  (48), 
equation  (53)  implies  a  starting  current  of 

I,  =  -  e(u2  ~1)/2  ,  u  <  0  .  (55) 

The  starting  current  is  shown  as  a  function  of  the 
frequency  variable  u  in  figure  5. 

Since 


equation  (54)  must  be  solved  self-consistently  for 
the  frequency  of  orotron  oscillation.  However,  if 
the  Q  of  the  open  resonator  is  sufficiently  large,  we 
may  replace  u  by 


on  the  right-hand  side  of  equation  (54).  Using 
equation  (55)  in  equation  ( 50)  and  substituting  the 
result  into  equation  (54),  we  obtain 

oj  =  (o0  X 


e^2 

(2tt)  i/J  Quo 


.  (56) 


The  second  term  in  equation  (56)  is  a  very  small 
correction  to  the  first,  since  Q  ~  5000  and 
Uy  ~  —  1 .  Therefore,  the  output  frequency  is  very 
close  to  the  resonant  frequency  of  the  passive  open 
resonator.  However,  this  correction  determines 
the  electronic  tuning  characteristics  of  the  orotron 
as  follows.  The  change  in  output  frequency  with  a 
change  in  accelerating  voltage  V  is 


do>  _  (d0  dv0  dc^,  d 

dV  _  (2tt)  1/2  Q  dV  dv0  di^ 
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Figure  5.  Variation  of  starting  current  in  orotron. 


we  obtain  (using  the  properties  of  the  confluent  Weconsidertheevaluationofequation(59)forthe 
hypergeometric  function)12  HDL  75-gigahertz  orotron  experiment.  Using 


doj 

dV 


W><0  0 

(8ir) 1/2  QVvoU2 


tu0  =  2tr  X  75  X  10g/second, 
wv  =  1  centimeter. 


We  are  primarily  interested  in  the  evaluation  of 
equation  (58)  in  the  region  near  Ug  =  —  1 ,  where 
the  starting  current  of  equation  (55)  is  a  minimum. 
The  hypergeometric  function  in  equation  (58)  has 
the  value  el/2  at  this  point;  and  so  equation  (58) 
becomes 


12  /  S.  Gradsfueyn  and!  M  Ryshik.  Table  of  Integrals.  Series,  and  Products. 
Academic  Press,  Sew  York  (1965). 


V  =  2500  volts, 

v0  =  3  X  107  meters/second, 

Q  =  5000  (estimated),  (60) 

we  obtain  for  the  electronic  tuning  transconductance 

df _ 1_  /  dtu  \  _ 

dV  “  2rr  \dV/0  “ 

0.31  megahertz/volt  (61) 


5.  CONCLUSIONS 

We  have  derived  a  linear  theory  of  electron 
bunching  in  orotrons  based  on  the  time  evolution  of 
the  distribution  function  representing  the  electron 
beam.  This  approach  to  bunching  has  the  ad¬ 
vantage  over  conventional  approaches6 1  that  it  is 
easily  generalizable  to  a  nonlinear  theory  in  which 
saturation  effects  are  included.  The  result  for  tht 
power  radiated,  which  is  derived  in  section  2, 
applies  not  only  to  the  orotron,  but  to  any  device  in 
which  a  longitudinal  electric  field  acts  as  the 
bunching  field.  Sections  3  and  4  consider  the 
evaluation  of  the  starting  current,  IJt  and  the 
electronic  tuning  transconductance,  df/dV.  These 
quantities  were  evaluated  for  the  planned  orotron 
experiment  at  75  gigahertz;  the  results  were 
Is  =  34.5  milliamperes  and  df/dV  =  0.31 
megahertz/volt.  It  is  remarkable  that  these  results 
agree  well  with  published  experimental  data  taken 
in  the  Soviet  Union.  (For  example,  in  one  experi- 


1  V  P  Sheslopalov,  Diffraction  Electronics.  Khar'  kov.  Moscon  i  I07fit 
(Trans,  by  IS.  Joint  Publications  Service.  April  1978). 

u  M  B  Txeytlin,  G.  A.  Bernashevskiy.  V  D  Kotov,  and  I.  T  Tsilson'.  Radio 
Eng  Elect.  Phys..  11(1977).  132. 

K.  Mizuno,  S.  Ono.  and  >’.  Shibata.  IEEE  Trans.  Electron  Dev..  FP-20 
(1973).  749. 


mental  device  similar  to  the  one  under  construction 
at  HDL,  values  of  Is  =  28  milliamperes  and 
df/dV  =  0.37  megahertz/volt  were  measured.’) 

An  extension  of  the  theory  to  the  nonlinear 
regime  is  necessary  to  calculate  the  actual  output 
power  of  the  orotron  as  a  function  of  the  electron 
beam  current.  The  formalism  introduced  in  this 
report  will  allow  such  a  calculation  to  be 
performed  in  a  straightforward  manner.  This 
nonlinear  theory  will  also  allow  us  to  consider  such 
effects  as  second  harmonic  generation  in  the 
orotron.  This  calculation  will  be  the  subject  of  a 
future  report  in  this  series. 

In  section  3,  we  consider  the  electromagnetic 
field  distribution  near  the  reflecting  diffraction 
grating  of  the  orotron.  In  particular,  the  unknown 
coefficients,  ar,  in  the  expansion  of  the  electromag¬ 
netic  field  are  determined  by  solving  a  boundary 
value  problem  near  the  grating.  We  have  per¬ 
formed  these  calculations  and  will  report  on  the 
results.  The  grating  calculation  allows  not  only  the 
determination  of  the  ar,  but  also  the  calculation  of 
the  quality  factor  Q  of  the  open  resonator  from  first 
principles. 
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